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SPECTRA OF Sol- MANIFOLDS: 
ARITHMETIC AND QUANTUM MONODROMY 

A.V. BOLSINOV, H.R. DULLIN, AND A.P. VESELOV 


Abstract. The spectral problem of three-dimensional manifolds admit¬ 
ting Sol -geometry in Thurston’s sense is investigated. Topologically are 
torus bundles over a circle with a unimodular hyperbolic gluing map A. The 
eigenfunctions of the corresponding Laplace-Beltrami operators are described 
in terms of the modified Mathieu functions. It is shown that the multiplicities 
of the eigenvalues are the same for generic values of the parameters in the 
metric and are directly related to the number of representations of an integer 
by a given indefinite binary quadratic form. As a result the spectral statistics 
is shown to disagree with the Berry-Tabor conjecture. The topological nature 
of the monodromy for both classical and quantum systems on Sol -manifolds 
is demonstrated. 


1. Introduction 

It has been known since the nineteenth century that in dimension two there is a 
close relationship between geometry and topology. Namely each compact orient able 
surface admits a metric of constant curvature: positive if it is a topological sphere, 
zero if it is a torus and negative if it has genus more than 1. 

In dimension three the situation is much more sophisticated. The major devel¬ 
opment here was due to Thurston m who put forward the famous Geometrisation 
Conjecture: any compact orientable 3-manifold can be cut by disjoint embedded 

2- spheres and tori into pieces, which after gluing 3-balls to all boundary spheres, 
admit one of 8 special geometric structures. These special 3-dimensional geome¬ 
tries are the standard Euclidean E 3 , spherical S 3 and hyperbolic H 3 geometries, 
the product geometries S 2 x R and H 2 x M and three geometries related to the Lie 
groups SL 2 (M.), Nil and Sol. 

The last group Sol is the 3-dimensional solvable Lie group, which is isomorphic 
to the group of isometries of Minkowski 2-space. The corresponding metric has the 
least symmetry of all the 8 geometries as the identity component of the stabiliser 
of a point is trivial. 

The structure of 3-manifolds admitting any of the seven geometries excluding 
the most complicated hyperbolic case H 3 is pretty well understood. In particular a 

3- manifold M possesses Sol -geometric structure if and only if M is finitely covered 
by a torus bundle over S 1 with hyperbolic gluing map. For all other 6 geometries 
M must be a Seifert fibre space (see e.g. pEO), so the SoLmanifolds are special 
from this point of view. 

Their special role in the theory of dynamical systems became clear after a recent 
paper [|| by Taimanov and one of the authors, who showed the surprising fact that 
although the geodesic flow on .SoLmanifolds is integrable in the sense of Liouville 
(but not in the analytic category) it has non-zero topological entropy ! 
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In the present paper we investigate the quantum version of the geodesic flow 
on SoZ-manifolds, which is the spectral problem for the corresponding Laplace- 
Beltrami operator A. We describe the spectra explicitly in terms of the spectrum 
of the modified Mathieu equation. These spectra are degenerate and have very 
interesting arithmetic. The multiplicities are directly related to the numbers of 
representations of a given integer by an indefinite binary quadratic form determined 
by the corresponding hyperbolic gluing map. This allows us to conclude that the 
spectral statistics for SoZ-manifolds is not Poisson contrary to the well-known Berry- 
Tabor conjecture. 

Note that the Sol -structure on SoZ-manifolds is not unique in the same way as 
the flat structure on a torus is. The spectra of tori are very sensitive to a change 
of the flat metric: if we change the periods slightly the degeneracy will essentially 
disappear. The fact that this does not happen with 5oZ-manifolds shows the rigidity 
of the spectra and can be considered as a reflection of the hyperbolicity hidden inside 
the topology of SoZ-manifolds. 

We should mention that a deep relation of S'oZ-manifolds with arithmetic was 
known before (see e.g. pun Ei). In particular, Hirzebruch m and Atiyah, Don¬ 
nelly and Singer P discovered a remarkable relation between topological ” signature 
defects” of SoZ-manifolds and arithmetical L-functions. 

From the dynamical point of view the arithmetic and topology reveal themselves 
through Hamiltonian monodromy [5j. Its quantum analogue - quantum monodromy 
- is a relatively new phenomenon PEI El, which still needs better understanding. 
An interesting feature of our case is that the corresponding grid of the quantum 
states can be described explicitly and nicely visualised (” Sol-flower”, see fig. EBP 
below). This is probably the first example of quantum monodromy of that kind. 

The structure of the paper is following. First we introduce the class of Sol- 
manifolds and describe the classical geodesic dynamics and the corresponding Hamil¬ 
tonian monodromy. Then we review the facts from classical number theory about 
the relations between binary quadratic forms and the modular group SL( 2, Z). In 
section 5 we consider the spectral problem for the corresponding Laplace-Beltrami 
operator and find the eigenfunctions in terms of the modified Mathieu functions. 
The arithmetic of the multiplicities of the eigenvalues is discussed in detail in sec¬ 
tion 6. The semiclassical analysis of the problem is done in section 7 in relation 
with Weyl’s law. In section 8 we discuss the spectral statistics in the context of 
the Berry-Tabor conjecture [2j. The quantum monodromy for SoZ-manifolds is 
discussed in the final section. 


2. SoZ-MANIFOLDS 

In this paper we restrict ourselves to the main class of S'oZ-manifolds, which 
are T 2 torus bundles over a circle S 1 with hyperbolic gluing maps with positive 
eigenvalues. More precisely, consider the action of Z on M 3 = T 2 x R generated 
by the following transformation Ta■ Let {x,y) be standard periodic coordinates on 
T 2 defined modulo 1 , and z £ (—oo, +00) be a coordinate on R. Then in these 
coordinates the transformation Ta is given by 

( anx + a 12 y 
a 2 \X + a 22 y 
z + 1 



(1) 
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where A = ^ ^ ^ J £ SL(2, Z) is an integer hyperbolic matrix, which defines 

a hyperbolic automorphism of the 2-torus. The corresponding SoZ-manifold M\ is 
defined as the quotient M 3 /Z by this action. 

Let A and A -1 be the eigenvalues of A and we assume that A > 1. The Sol- 
manifolds with negative A are covered by those with positive eigenvalues. 

Together with ( x,y,z ) we shall use another coordinate system (u,v,z) on . 
where (u, v) are linear coordinates on the fibres related to a positively oriented 
eigenbasis of A. The transformation Ta in these coordinates is given by 

Xu \ 

\~ 1 v 

z + 1 J 



One should note that unlike (x,y), the new coordinates (u,v) are not periodic 
on the tori T 2 anymore: two pairs (u, v), (v !, v') define the same point on T 2 if and 
only if (u — u',v — v') = k(c\,c\) + m(c 2 ,c|), where k,m £ Z and e\ = ( c\,c 2 ), 
e 2 = (cl, cl) is the basis of the lattice T associated to T 2 : 


{ an ai2 

V °2i a22 



The Rienrannian metrics on S'oLmanifolds come from right-invariant metrics on 
the universal covering of M\, which has the natural structure of a solvable Lie 
group Sol. Topologically this group is R 3 with a multiplication of the form 


(u, v, w) * ( u', v', w 1 ) = (u + e w v!, v + e w v ', w + w'). 


One can realise it as the group of 3 x 3 matrices of the form 

/ e w 0 u 
0 e~ w v 

\ 0 0 1 

The SoLmanifolds M\ we consider are the quotients of the group Sol by the 
discrete subgroups Ga corresponding to w = min A, m € Z and ( u,v ) = ke i + le 2 
belonging to the integer lattice T described above, z = w/ In A. 

The right-invariant metrics on the group Sol correspond to the following class 
of metrics on the SoZ-nranifold Ma- 



( 3 ) 

where 


ds 2 = a(z) dx 2 + 2f3(z)dxdy + r )(z)dy 2 + dz 2 


a(z) P(z) 
f3(z) 7 (z) 


= exp(— zB) 1 


a (3 

(3 7 


exp(— zB). 


Here a, (3 ,7 are real parameters with the only condition that the form ds 2 = adx 2 - 
2f3dxdy + yd y 2 is positive definite and B is defined by the relation exp B = A: 


B = 


-1 


In A 
0 


0 

— In A 


One can consider a more general metric allowing a constant coefficient at dz 2 but 
this will lead only to a general scaling. 
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3. Geodesic flows on SoZ-manifolds: integrals and Hamiltonian 

MONODROMY 


Thus, the Hamiltonian of the geodesic flow on M\ in (u, v, ^-coordinates can 
be written as 

H = I (Ee 2zlnX p 2 u + 2F PuPv + Ge~ 2zlnX p 2 ) + ±p 2 z , 

where E, F, G are real parameters: E > 0, G > 0, EG — F 2 > 0. It is invariant 
under the following transformation 


/ u \ 


/ Xu \ 

V 


A -1 v 

z 


z+1 

Pu 


A ~ 1 p u 

Pv 


A p v 

\ Pz / 


\ Pz ) 


and, of course, under the translations by the elements of the lattice T. The same 
property must be satisfied for any smooth function on T* in particular, for the 
first integrals of the geodesic flow. 

Since FI depends neither on u, nor on v, the corresponding momenta p u and 
p v are local first integrals of the geodesic flow. However, being not invariant un¬ 
der 0 , they are not well defined on the cotangent bundle T*M\. That is why, 
to get global first integrals, we need to replace p u , p v by two smooth functions 
h{Pu,Pv), f-2(Pu,Pv) invariant under the transformation ( p u ,p v ) —> (A _1 p M , Ap„) 
(or, speaking in more general terms, by the invariants of the Z-action on the cotan¬ 
gent plane generated by the hyperbolic linear transformation A T 1 ). 

One invariant function is evident: Q = p u p v . To find another one we introduce 
the following expression which will be useful also in the future 


( 5 ) 


a = 



2 In A 


Under the transformation 0 a changes in a very simple way: 


ot{p u ,Pv) —► a(p u ,Pv) - 1 


Thus, as a second integral we can take any function of a with period 1, for 
instance, cos(27ra) or sin(27ra). However these functions are not smooth at p u = 
p v = 0. To avoid this difficulty and to get the first integrals in a more symmetric 
form we put: 

/i = R(Q) cos27ra, 
fi = R{Q) sin27ra, 

where 

R(Q) = vWI ex p(-^2)- 

Remark. The fact that the second integral is not analytic is not accidental: 
the theorem proved by Taimanov m implies that ^oZ-manifolds do not admit 
integrable geodesic flows with analytic integrals (see |J] for more details). 

We are going to show now that one can see the topological structure of the Sol- 
manifolds by looking at the Hamiltonian monodromy of the geodesic flow. For that 









SPECTRA OF So/-MANIFOLDS 


5 


we will have to investigate the bifurcation diagram (i.e. the set of critical values) of 
the momentum mapping restricted to the isoenergy surface E\ = {H = 1}: 

(6) E A = (f 1 J 2 ):E 5 A ->R 2 . 

Proposition 1. The bifurcation diagram of the momentum mapping T A consists 
of two circles 

fl + fl = R\Ql), 

where Q*± = (F ± \fEG)~ 1 , and the point (0,0), the centre of these circles. The 
set of critical points consists of five parts: a) four one-parameter families Li (i = 
1,... ,4) of (degenerate) 2-dimensional tori lying in the cotangent bundle and given 
by (a is a parameter): 

z = —a, u and v are arbitrary, 


0 2ot In A 


Pz = 0, Pu = ± 


W 


i + 


Pv = ± 


and 


Veg 

z = —a, u and v are arbitrary, 


0 2a In A 


Pz = 0, Pu = ± 


W 1 - 


, 


Pv — T 


g— 2a. In A 

G 

(l , F ) 

Veg) 

g— 2a. In A 


Veg) 

b) the critical set N given by the equation Q = p u p v = 0. 


\ G l 1 VTg) 


Proof. We are interested in the singularities of T A or, which is the same, those 
of the Liouville foliation. These singularities can be of two types. To explain their 
nature we first consider the geodesic flow on the covering manifold M . On this 
(non-compact) manifold the integrals of the flow are simply p u and p v . Consider 
the Liouville foliation for this covering system. Its singular leaves correspond to 
the critical points of the momentum mapping 

F = (Pu,Pv) ■ E 5 -> K 2 , 

where E 5 = {E[ = 1} C T*M. Obviously, these leaves remain singular after the 
natural projection E 5 —» E\. These are singularities of the first type. 

On the other hand some new singularities appear since instead of p u and p v we 
have to consider more complicated functions fi and f 2 . In other words, these are 
singularities of the map (p u ,p v ) (/i, / 2 )- 

Let us treat both cases in turns. It is easily seen that p u and p v are function¬ 
ally dependent, as functions on E 5 = {H = 1} if and only if two conditions are 
simultaneously satisfied: 1) = 2 p z = 0 and 2) ^ = 0. Taking into account the 

condition H = 1, we obtain a system of equations 

Ee 2\o g Xz p 2 u _ Ge -21ogAzp2 = 0; 

Ee 2 \o g \z p 2 u + 2 F PuPv + Ge _21og Xz pl = 2 
The first equation gives 
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Figure 1. The topological structure of the singular leaf is M\ x K 
where K is given in the figure. 


Now solving this system with respect p u and p v (after substituting 2 
find four distinct solutions: 


1) Pu 

2) p u 

3) p u 

4) Pu 



Pv 

Pv 

Pv 

Pv 



—a ), we 


Thus, for each value of a we obtain four 2-dimensional invariant tori in T*M\. 
All of them are diffeomorphically projected onto the same T 2 -fibre T 2 a = {z = 
const = — a} C M. Varying a, we obtain 4 families of degenerate Liouville 2-tori 
Li, i= 1,... ,4. 

It is easy to verify that for each family Li the value of Q = p u Pv is constant and 
equal to Q+ = (F + VEG ) _1 for L\ and L 2 , and equal to Q*_ = (F — \/£G ) _1 
for L 3 and L 4 . Hence the image of L 1 and L 2 is the circle f 2 + f 2 = R 2 (Q* + ), 
and analogously the image of L 3 and L 4 is the other circle f? + f 2 = R 2 (Q*_), as 
required. 

The singularities of the second type come from those of the mapping (p u ,Pv) —► 
(/l,/^)- It can be easily seen that the critical points of this mapping are defined 
by the equation Q = p u p v = 0. This implies immediately /1 = fi = 0 which gives 
a single point on the bifurcation diagram, namely the centre of the circles. 

Notice that topologically the subset N = {Q = p u p v = 0} C T*M\ is homeo- 
morphic to the direct product M\ x K, where K is a graph that consists of two 
vertices and four segments connecting them (see fig. This follows immediately 
from the parallelizability of M\ and the simple observation that in each cotangent 
space the conditions p u Pv = 0, H = 1 define a graph homeomorphic to K. 

Now we are able to describe the global structure of the foliation of the isoenergy 
surface E\ into Liouville tori. 

If we remove the singular set from the isoenergy surface we obtain four families 
of 3-dimensional Liouville tori distinguished from each other by signs of p u and p v : 

a) p u > 0 ,p v > 0 ; 

b) p u < 0 ,p v < 0; 

c) p u > 0 ,p v < 0 ; 

d) p u < 0 ,p v > 0 . 
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The families a) and b) are isomorphic (more precisely, they transform into each 
other by the globally defined time reversal automorphism of the geodesic flow 
(u,v,z,p u ,p v ,p z ) -> {u,v,z,-p u ,-p v ,-p z )). 

The same is true for the families c) and d). 

Each Liouville 3-torus is uniquely determined by the values of two integrals Q 
and a mod 1, where the values of Q form the interval (0, Q+) in the first two cases 
and (QI,0) for the other two cases. In particular, in each of the cases, the base 
of the T 3 -foliation is homeomorphic to a punctured disc. As Q —> 0, the Liouville 
torus approaches the singular set N. As <3 —> Q±, the torus shrinks into one of the 
degenerate 2-tori described above. 

Thus, the base of the global Liouville foliation on E\ = {H = 1} can be consid¬ 
ered as four discs glued together at their centres. All interior points of these discs 
except the centre correspond one-to-one to regular 3-dimensional Liouville tori, the 
boundary circles of the discs correspond to the families L, of degenerate 2-tori, and 
finally, the common center of the discs corresponds to the singular set N. □ 

The image of each family under the momentum map is a 2-disc with the cen¬ 
ter removed. This is exactly the situation when we can talk about Hamiltonian 
monodromy 0. 


Theorem 1. For each family of Liouville 3-tori there exist a basis of cycles in 
the first homology group of the tori in which the Hamiltonian monodromy has the 
matrix 

A 0 
0 1 


Proof. This fact can be observed in many different ways. We shall follows the 
definition of Hamiltonian monodromy and will explicitely compute the deformation 
of Liouville tori and the final gluing map. 

Consider an arbitrary Liouville 3-torus T 3 = Tq 0 ao - In coordinates, this torus 
is given by three conditions: 


(7) 


Ee 2z ln x p 2 u + 2 Fp u p v 
Q(p u ,Pv) — PuPv — Qo 

Ot{p u ,Pv) = ^ 


Ge -2z\n\ p 2 


= 2 , 


= Oq mod 1 


More precisely, these conditions define a disjoint union of two or four tori, which 
differ from each other by the signs of the momenta p u and p v . We consider one of 
them Tq q ao by putting for definiteness p u > 0, p v > 0. 

For our purposes first we need to explain why the above conditions define indeed 
a three-dimensional torus and to describe the basic cycles on this torus. Notice that 
the common level set 0 of the first integrals can be regarded from two slightly 
different points of view: as a subset in T*M 3 and that in T*M\. However one can 
show that the natural projection T*M 3 —> T*M\ restricted to this level set is a 
diffeomorphism (no points are glued between them). Thus, in fact there is no real 
difference between these two points of view. In particular, instead of conditions 
PuPv = Qo, oi(p u ,p v ) = a 0 mod 1 we may simply assume that the momenta p u , p v 
themselves are constant. Then the conditions 0 can be rewritten as: 


p u = const , p v = const , u and v are arbitrary, 


and 

( 8 ) 


ci cosh(2 ln X(z + ccq )) + p 2 z = C 2 
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where Ci = 2\/EG\Qo\, C 2 = 2 — 2 FQ 0 . We see that the variables separate and the 
fact that this systems defines a 3-torus becomes evident. Indeed, the variables u,v 
’’run” over a two-dimensional torus and the last equation defines a simple closed 
curve on the plane M. 2 (z,p z ). In other words, we have a natural splitting of Tq q ao 
into the direct product T 2 x S 1 . Thus, as basic cycles on Tq we can take the 
cycles on T 2 (u,v) related to the original coordinate system (x,y) (see above) and 
the third cycle defined by ©■ 

Now let us look at what happens to this torus if we change the parameters Qo 
and «o in such a way that the point Ta{Tq 0 ao ) moves inside the image of the 
momentum mapping around the singular point !Fa{N) = (0,0). It is easy to see 
that this deformation just means that we change the value of a , while Q can be 
chosen to remain constant: 


Q(t) = Q o, a(t)=ao+t, t G [0,1]. 


Consider the family of mappings 

(l>t(u,v,z,p u ,p v ,p z ) = (u, v, z- t, e tlnX p u ,e~ tlnX p Vl p z ). 


It is not hard to see that the image of Tq o under </> t is exactly Tq q ao+t 
and (j>t : Tq o ao —> Tq o ao+t is a difeomorphism. In other words, <j>t defines the 
deformation of Liouville tori we need. 

At the moment t = 1 the torus comes back to the initial position, i.e., Tq q ao = 
Tq 0 a()+ i, and we obtain the monodromy map 

: T Qo,<xo “* T Qo,oc 0 = T Q 0 ,a 0 +l 


Now our goal is to describe the corresponding automorphism of the first homol¬ 
ogy group: 

tu'-H 1 (T^ ao )=Z 3 ^H 1 (T^ ao )=Z 3 . 

Using the identification 0 we see that the map <p t can be rewritten as follows: 


< U ^ 


/ Ait \ 

V 


A” 1 ?; 

z 


z 

Pu 


Pu 

Pv 


Pv 

\Pz ) 


V Pz J 


We see that the only transformation is related to the variables u and v. Moreover, 
this transformation is exactly the original hyperbolic automorphism A : T 2 —> T 2 . 
Taking into account the natural splitting Tq 0 ao = T 2 (u,v ) x S 1 (z,p ~) we conclude 
immediately that the monodromy matrix in the chosen basis is 


A 

0 



□ 

We conclude this section with the discussion of the geodesics on SoZ-manifolds. 
They have different properties depending on the types of leaves of the Liouville 
foliation which they belong to. 

First consider the geodesics lying on Liouville tori of dimension three. They 
are characterized by the property that all momenta p u ,Pv and p z differ from zero. 
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More precisely, the signs p u and p v always remain the same, whereas the sign of p z 
changes. This happens when 2 reaches the value 


z± 


±cosh 1 


( Fp u p v \ 
\VEG\p uPv \) 


2 In A 


a{p u ,Pv)- 


Two levels z = z+ and z = Z- are exactly the caustics of the Liouville tori that 
contains a given geodesic. The situation is quite similar to that on a surface of 
revolution where the motion takes place between two levels of 2. 

It is easy to see that the distance between these levels z + — tends to infinity 
as p u p v tends to zero. From this it follows that the corresponding geodesics rotate 
many times (along the base S' 1 ), then turn back, after this go in the opposite 
direction, then turn back and so on. As p u p v tends to zero the number of rotations 
in one direction until turning back (or, which is the same, between two caustics) 
increases up to infinity. 

If PuPv = 0, then we are on the singular level. The corresponding geodesics have 
the following behaviour. If both p u and p v vanish, then we obtain the family of 
geodesics 

u = const , v = const , z = t. 


Such geodesics obviously form an invariant submanifold N + in T*M which is 
diffeomorphic to M. Exactly on this submanifold the geodesic flow is chaotic and 
has positive entropy. Indeed, the time-one map transform each fibre T® into itself 
by means of the hyperbolic automorphism A. As well known, the entropy of A : 
T 2 -> T 2 is In A > 0. 

There is another invariant submanifold _/V_ with the same properties formed by 
vertical geodesics going in the opposite direction: 


u = const , v = const, z = —t. 

From the viewpoint of the ambient geodesic flow N + and N— are hyperbolic 
invariant subsets. The stable manifold coresponding to N + is given by p v = 0, the 
unstable one is p u = 0. For N_ the stable and unstable manifolds interchange. 
The geodesics satisfying the condition p v = 0 as t —> +00 asymptocally approaches 
N+, in particular, p z — > +1. But there is t = to when p z changes sign so that for 
t —> —00 the geodesic approaches to IV_. The geodesics satisfying p u = 0 behave 
in the opposite way. 

In slightly other terms this structure can be described as follows: there are 
two hyperbolic submanifolds diffeomorphic to M\, they are connected by 4 four¬ 
dimensional separatrices, see fig. ^ 

Finally we would like to mention an interesting phenomenon which one would 
not expect from an integrable geodesic flow on a compact manifold. Namely, one 
of the action integrals diverges as the integral Q —> 0 with the energy fixed (see 
the calculations and footnote in Section 7). Normal scenario would be when ap¬ 
proaching the singular level some of the cycles of the Liouville tori shrink so the 
actions will stay finite. The fact that this not true for .SoZ-manifolds when one ap¬ 
proaches the singular (chaotic) level demonstrates once again the peculiar nature 
of this system. 

To discuss the quantum case we will need some facts from the classical number 
theory, which we present in the next section. 
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4. SL( 2,Z) AND BINARY QUADRATIC FORMS 

The content of this section is well-known (see e.g. HHUmED). 

Let A = ( 0,11 0,12 ] £ SL( 2,Z) be an integer hyperbolic matrix. Hyperbol¬ 

ic 021 «22 J 

icity as before means that its eigenvalues are real and distinct. We would like to 
consider A as the automorphism of the lattice C = Z © Z £ R 2 by choosing some 
basis e \ 1 e 2 in this lattice. 

For any such A we can define the following integer binary quadratic form Qa by 
the formula 


(9) Av A v = Qa(v)ci A e 2 , 
where v is a vector from R 2 . Explicitly if v = xe\ + ye 2 then 

(10) Qa(x, y ) = det ( ° llX + ai2V X \ = -a 21 x 2 + (an - a 22 )xy + a 12 y 2 . 

\ &2ix + a 22 y y J 

It is easy to see from the definition that this form is invariant under the action of 
A: 

Qa{Av) = Qa(v). 

Notice that Qa has the discriminant 

D = (an — a 22 ) 2 + 4ai 2 a 2 i = (an + a 22 ) 2 — 4(ana 22 — ai 2 a 2 i) = (an + a 22 ) 2 — 4, 
which is exactly the discriminant of the characteristic equation of A: 

A 2 — (an + a 22 )A + 1 = 0. 

In particular, since A is hyperbolic the form Qa is indefinite. Note that the dis¬ 
criminant D cannot be a total square. 

In general the coefficients of the quadratic form Qa may have a common factor. 
Let 

(11) Q A (x,y) = ax 2 +bxy + cy 2 


be its primitive form after division of Qa by the largest common factor. It is defined 
correctly only up to a sign. 

Thus to each integer unimodular hyperbolic matrix A we relate an indefinite 
integer primitive quadratic form Qa- 

Conversely, suppose we have such a form Q(x, y) = ax 2 + bxy + cy 2 . We would 
like to describe all A from SL( 2,Z) which preserve this form. Such A are called 
the automorphs of Q. Let 

d = b 2 — 4ac 


be the discriminant of Q which we assume not to be a total square and consider 
the corresponding Diophantine equation called Pell's equation: 

(12) X 2 - dY 2 = 4. 


Then the group of automorphs consists of matrices of the form 


A = ± 


X-bY 

2 

aY 


-cY 

X+bY 


2 
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where ( X, Y) are the solutions of Pell’s equation. Modulo ±J this group is cyclic 
with generator 


(13) 


A 0 = 


—cYq \ 

aY 0 Xo±bYo ) > 


where (Xo,lo) is the fundamental solution of this equation. 

Recall that (Xo, Yq) is the fundamental solution of Pell’s equation if Xq > 0, lo > 
0 and Xo + VdYo is minimal among all such solutions. The classical result about 
Pell’s equation says that all other solutions can be found from the relation 

X + VdY / X 0 + VdY 0 

2 l 2 



where n = 0,1,.... One can find the fundamental solution from the continued 
fraction of \fd. This structure of the solutions of Pell’s equations induces the cyclic 
group structure for the automorphs. 

Notice that the form Qa corresponding to the matrix m has the form 

Q = Y 0 (ax 2 + bxy + cy 2 ). 

Let us call a hyperbolic element A from SL( 2,Z) primitive if it can not be 
represented as a power of any other element from SL( 2,Z). 

Thus we have described a natural correspondence between the primitive binary 
indefinite forms Q and primitive elements A from SL( 2,Z). In particular, it helps 
us to answer the question if a given integer unimodular matrix A is a primitive or 
if not which power of a primitive matrix it is. 


5. Spectrum and eigenfunctions of the Laplace-Beltrami operator 
Let us now discuss the quantum geodesic problem on the SoZ-manifold M\: 

(14) —A'if = £if, 

where A is the Laplace-Beltrami operator on M\ and if = il)(P,£),P £ M\. In 
coordinates (u, v, z ) the Laplace-Beltrami operator has the following explicit form: 


(15) 


A = Ee 2zlnX 


d 2 d 2 

du 2 ^ dudv 


+ Ge~ 2zlnX 


d 2 d 2 
fo 2 + dz 2 ' 


This is a self-adjoint operator in the Hilbert space L 2 {M\) where the integration 
measure on m a is induced by the Riemannian metric ©• In both (x,y, z) and 
{u, v, z) coordinate systems the corresponding measure dy is proportional to the 
standard Lebesgue measure on R 3 . 

Because the coefficients of A depends only on z it is quite natural to separate 
variables and look for the eigenfunctions of A of the form 




where 7 is an element of the dual lattice T* corresponding to the T 2 -fibres and 
w = (u,v) (so the scalar product is defined modulo Z). 

By substituting into the Schrodinger equation o, tUt we get 


(S“ 87r2 ^ |Q(7)l ( cosh ( 2 ln Kz + <*(7))) + FSg }§^ ] ) /) e27ri(7,W) > 
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where Q{ 7) = ('y,e u )( / y,e v ) is a quadratic form on the lattice F*, and 



Here e u and e v are the eigenvectors of A related to the eigenvalues A and A -1 
respectively and the basis e u , e v is assumed to be positively oriented. Notice that 
a is the same as before in 0 if we replace p u by (7, e„) and p v by (7, e v ). 

To clarify the meaning of the coefficient in front of the cosh let us consider the 
basis e*,e* in R 2 * dual to e u ,e v . The vectors e* and e* are also the eigenvectors 
of A* with the eigenvalues A and A -1 respectively. By definition we have 7 = 
(7,e u )e* + (7,e„)e*. Since <5(7) is obviously invariant under the action of A* it is 
natural to compare it with the binary form Qa * defined in the section 3 . We have 

- 4*7 A 7 = (A( 7 , e u )e* + A _1 (7, e v )e*) A ((7,e u )e* + (7,e„)e*) 

= ( 7 , e„)(7, e„)(A - A" 1 ^ A e*. 


Let lii I2 be a positively oriented basis in the dual lattice F*, then by definition 
Al*7 A 7 = Q a * (7)/iA! 2 . 

From these calculations and from the equalities E = |e*| 2 ,G = |e*| 2 it follows 
that 

yfEG\Q{l)\=c\Q A *{ 1 )\, 

where 


( 16 ) 


c = c{A- E, F, G) 


A{D*) 
\J~D sin 9 


1 

VDA(T 2 ) sin 8 


^(□*) is the area of the dual basic parallelogram n(e*, e^) (which is the inverse of 
the area of the fibre T 2 ), D = (A — A -1 ) 2 is the discriminant of the characteristic 
equation of the matrix A (or equivalently A*), and 9 is the angle between e* and 
e*. Thus we have proved the following 


Proposition 2. A function T = e 2m ^’ y,w ^ f(z) satisfies equation \ 14 \) if and only if 
f(z) satisfies the modified Mathieu equation 

( 17 ) + W(l)\ cosh 2 p(z + 0(7)^ f(z) = A f(z), 

where p = In A, 1/(7) = 87 t 2 cQa*( 7 ) and 0(7) is given above. The eigenvalues £ 
and E are related by the shift 

( 18 ) £ = A + 1/(7) cos 6 . 


Recall that the modified Mathieu equation is the cosh-version of the standard 
Mathieu equation 

d 2 y 

—- + (a cos 2 fix + b)y = 0. 
ax 2 

Its solutions are known as modified Mathieu functions (see e.g. .23 ECO)- They 
appear also in the theory of Coulomb spheroidal functions da, where one can find 
some related numerical results (see also USD- 

Let A = Afc(i/), k = 1 , 2 ,... be the spectrum of the corresponding modified 
Mathieu operator 

d 2 

M = + \v\ cosh2 pLZ 

d z 2. 








SPECTRA OF SoJ-MANIFOLDS 


13 


and fj,k(z) be the corresponding solutions of 

Thus, to each element 7 of the dual lattice F* we associate the functions ’F 7 i a,(z(, v, z 
e 27 ri ( 7 ’“')/ 7 i fe(^). The problem with these functions is that they are well defined on 
the covering space M 3 = T 2 xI but not on the SW-manifold M\ itself because 
they are not invariant with respect to the transformation £EJ, 0. One can try to 
construct the genuine eigenfunctions of A on by averaging these functions with 
respect to the action of Z on M 3 generated by this transformation. It turns out 
that the averaging procedure works. 

To show this let us consider instead of \l/ 7 i fc(it, v, z) the following sum 

(19) $ 7ifc = ^2 X~ n v, z + n) =^2 ^A*" 7 ,fc(w, v, z). 

n£Z n£Z 

Because of the fast decay of the eigenfunctions / 7 ,fe(^) this sum is absolutely con¬ 
vergent. It is easy to see that it defines a well-defined function on M\ , which is an 
eigenfunction of the Laplace-Beltrami operator A. 

The eigenfunctions d> 7 i / c (u, v, z) on M\ actually depend only on the orbits [7] = 
{A* n (p/)} n( z% with respect to the action of A* on T*: $ lt k(u,v,z) = $[ 7 ] ifc (u, v, z). 

We should also consider separately the eigenfunctions related to 7 = 0. It is easy 
to see that the corresponding eigenfunctions have the very simple form 

(20) $o,s = 1, cos 27 tz, sin27r,z, cos 47 rz, sin 4772,... , cos2fc7T2, sin2fc7rz,... 
with the eigenvalues £k = (27r) 2 fc 2 . 

Theorem 2. The eigenfunctions of the Laplace-Beltrami operator < Fr 7 i *,( 16 , v, z), [ 7 ] £ 
T* \ {0}/A* and 4>o, s (-z) form a complete basis in Lz^M^). 

Proof. The independence and orthogonality of these functions are obvious. The 
only thing we have to verify is the completeness. To prove this we need to show 
that any smooth function 4> : M\ —> R which is orthogonal to each eigenfunction 
from the list is, in fact, zero. 

Consider such a function $(ui, 2 ) on the covering space M 3 and expand it as a 
Fourier series (with respect to w): 

$(w,z) = e 2 ”^’ w) a 7 {z) 

7 g r* 

with some smooth coefficients a 7 (z),z £ M. 

Lemma 1. For all 7 / 0 the functions a- t {z) have fast decay at infinity and thus 
belong to L 2 (R). 

Proof Since 4> is invariant with respect to the transformation 0 , we have 
<I>(z< 7 , z) = *I>(Hu>, z + 1). Hence 

e 2TTi( 7 ’ w )a 7 (z) = J2 e 27ri ^’ Aw) a 1 ( Z + 1) = ^ e 2 ^ A ^’ w) ai {z + 1). 

76 r* 7 gr* 7 er» 

Thus the Fourier coefficients satisfy the following property: 

a 7 (z T 1) — u^* 7 (z), 

or, more generally, 

a 1 {z + n) = aA* n -y(z ), n £ Z. 
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Since the Fourier coefficients a 7 of a smooth function decay fast for large 7 and 
A* k "/ for 7 7^ 0 tends to infinity we see that the functions a 1 {z) decay very fast and 
thus belong to L 2 (R). □ 

Now suppose that &(w,z) is orthogonal to the eigenfunction 3>[ 7o ],fc(u, v, z) = 
Y^ 7 i£Z®A* n 'y 0 ,k(u,v,z). Since the measure on M\ is proportional to the standard 
Lebesgue measure dudvdz we have 


0 = ($(w, z), $[ 7o ] lfc (>, z)) = / $(w, z)$ [7o]jfc (w, z)da 

= f 1 I E E / e 2 ^^e~ 2 ^ n ^dudv ) a y (z)f A *n JOtk (z)dz 

y 7 er* n& Jt2 ) 

= / (V / e 2 ^ n ^e~ 2 ^ A,n ^dudv\a A *n 10 (z)f A *n^ tk (z)dz 

J o \^z J t2 J 

= A{T 2 ) f Y a A -"-'y 0 (z)f A >n'y 0 ,k{z)dz . 

Jo 


We now use the property that f-y,k(z+n) = /A* n 7 ,fc(~) and a 7 ( z+n ) = a A * n j{z), 
Z to conclude that 



(z)f A ^ JO ,k(z)dz = 


n G 


/•1 _ p +00 

/ °To (z + n)f 70tk {z + n)dz= / a l0 (z)f l0tk {z)dz. 
"'° n£Z J -00 


Thus, the Fourier coefficients a 7o (z) for 70 7^ 0 belong to L 2 (R) and at the same 
time are orthogonal to all the functions / 7o ,fc(-z) which form a complete basis in 
L 2 (M). Hence for 70 7^ 0 the coefficients a l0 {z) = 0. 

This means that the function <f> must be of the form <F(w;, z) = a(z), where a(z) is 
periodic with period 1. Now using orthogonality to the functions m we conclude 
that a(z ) must be identically zero. □ 


Corollary 1. The spectrum of the Laplace-Beltrami operator on Sol-manifolds 
consists of two parts: the trivial part 

£ = £ k = 4/c 2 7t 2 , k = 0,1,... 

corresponding to the eigenfunctions Iffll and the non-trivial part 

£ = £z,[ 7 ] = MKM)) +K[ 7])cos6>, 1 = 1,2..., [7] G F* { 0 }/ A* 

related to the modified Mathieu equation m- 

The multiplicities of the trivial eigenvalues are 2 except for the ground state 
£ = 0 which has multiplicity 1. The multiplicities of the non-trivial part of the 
spectrum are much more interesting and the answer depends on the arithmetical 
properties of the gluing map A. We discuss this in the next section. 


6. Multiplicities of the eigenvalues and number theory 

As one can see from the previous section the eigenvalue of $[ 7 ] k (u, v, z ) depends 
on 7 only via Q A *( 7)- Thus the calculation of the multiplicity (for generic values 
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Py 



Figure 2. Fundamental domain of the lattice in ( p x ,p y ) with 
\Q\ < 30 2 for the cat-map. The hyperbola Q = — ll 2 illustrates 
the first example of a non-trivial degeneracy. 


of the parameters to avoid additional accidental coincidences) is reduced to the 
classical number theoretic problem of finding the number Nq(h) of integer solutions 
of the equation Q{x , y) = ax 2 + bxy + cy 2 = n for a primitive indefinite quadratic 
form Q different modulo its automorphs. Figure 3 illustrates this for the cat-map 
A with Q = — x 2 + xy + y 2 . 

For forms Q with certain discriminants there exists an effective formula which 
allows to compute Nq(u) To be more precise we need the following notion. We 
say that two forms Q and Q' are equivalent if there exists a transformation from 
SL( 2,Z) mapping one into another. It is easy to see that two equivalent forms 
must have the same discriminant d = b 2 — 4 ac. The converse is not true: there can 
be more than one non-equivalent forms with the same discriminant. 

Let h(d) be the number of classes of primitive forms with discriminant d. Note 
that the discriminant d = b 2 — 4 ac is always 0 or 1 modulo 4 and we assume as 
usual that it is not a total square. 

Remark. One should distinguish h{d) and the class number of ideals in the 
quadratic number field Q(\fd). They coincide only if the so-called negative Pell 
equation 


X 2 - dY 2 = -4 
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has a solution; otherwise h(d) is twice as big (see e.g. [Tl] . Chapter 16). The 
last property can be reformulated in terms of the period of the continued fraction 
expansion of y/d, but a more explicit description is unknown. 

If h(d) = 1 then all forms with the discriminant d are equivalent. In that case 
there is the following remarkable formula for the number Nq (n) when n is positive 
and coprime with d: 

(21) N Q (n) = N d (n) = E (f) ’ 

k\n ' ' 

where the sum is taken over all divisors of n and (^) is the standard Kronecker 
symbol (see Landau mi, Chapter IV.4). The Kronecker symbol is a real character 
modulo d , which has the following properties determining it uniquely: 

(1) If d and k are not coprime then (^) =0; 

(2) If d and k are coprime then (|) = ±1; 

(3) (*) = (f) (f) 

(4) for p odd prime which is not a divisor of d coincides with the Legendre 
symbol , which is 1 if d is quadratic residue modulo p and —1 otherwise; 

(5) (|) is 1 if d has residue 1 modulo 8 and —1 if it has residue 5 modulo 8. 

For its computation one can use the celebrated Law of Quadratic Reciprocity: if 
p, q are coprime positive odd numbers then 



Here is the list of the discriminants d up to 100 with h{d) = 1, see m 

5,8,13,17,20,29,37,41,52,53,61,65,68, 73,85,89,97. 

It is believed that there are infinitely many fundamental discriminants with h(d ) = 
1, but it is still an open problem. Notice that for positive definite forms it is known 
that there are only 9 fundamental discriminants with h(d) = 1 as it was conjectured 
by Gauss, namely 

d = -3, -4, -7, -8, -11, -19, -43, -67, -163. 

In general if h(d) > 1 the right-hand side of the formula 12111 gives the total 
number of representations of n by all non-equivalent forms with discriminant d. An 
interesting case is when the ideal class number of d is 1 but h(d ) = 2. In that case 
we have only two non-equivalent forms with discriminant d: Q and — Q and the 
formula <ED gives the number of the solutions of the equation |Q| = n. The first 
corresponding discriminants are: 

12,21,24,28,32,33,44,45,48,56,57,69,72,76,77,80,84,88,92,93 

(see HD- The only discriminants < 100 not listed in either table above are 
40,60,85,96 with h(d) = 2,4, 2,4, respectively. Note that most of these discrim¬ 
inants d are not of the form D = t 2 — 4, but they can still be obtained from 
A € 6X2(21) because D/d may be an arbitrary square. 

Now we are ready to describe the multiplicities of the eigenvalues £i,\y\- First we 
should take into account that the gluing map A £ SL(2,Z) and the corresponding 
form Q = Qa * may be non-primitive. Let us define the positive integers r = r(A) 
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and l = 1(A) from the relations A = Aq and Qa * = IQ, where Aq £ SL( 2,Z) and 
Q = Qa * are primitive. 


Theorem 3. The multiplicity m of the eigenvalue of the Laplace-Beltrami 

operator A for generic values of the parameters in the metric m is 

m( 7) = 2 r(A)N Q ,(n), 

where n = Q*( 7) = 1(A)~ 1 Qa* (7). When the discriminant d of the form Q* has 
class number 1 and n coprime d then Nq * ( n) can be computed using the formula 

m- 


Examples. 

1. Let 



1 

1 


be the so-called cat-map. Then A* = A and Qa = Qa * = — (x 2 — xy — y 2 ) are 
both primitive. The discriminant D = d = 5 has class number 1, so one can use 
the formula m to compute the multiplicity of the corresponding Sim- One can 
check that this leads to the formula 


to = 2(iV ±1 (n) - N ±2 (n)), 


where N±i(n) and N± 2 (n) are the numbers of divisors of n = Qa( 7) which have 
respectively the residues ±1 and ±2 modulo 5. 

This example shows that the multiplicities of the eigenvalues can be as big as 
we like: for example for n = 11 M the multiplicity is M + 1; for n a product of M 
distinct primes all ±1 mod 5 the multiplicity is 2 M . 

2. The matrix 


A = 


1 3 
3 10 


corresponds to d = 13 with D = 3 2 d and Qa * = —3(a; 2 — 3 xy — y 2 ) so 1(A) = 3. 

3. For 

' 5 2 


A = 


2 1 


we have d = 8 with D = 2 2 d and Qa * = Qa = —2(a; 2 — 2 xy — y 2 ) so 1(A) = 2. 
This example shows that 1(A) in general is not directly related to the largest square 
divisor of D. 

4. For 

' 1 3 

1 4 


A = 


d = D = 21 and Qa * = ~(3a; 2 + 3 xy — y 2 ), 1(A) = 1. Here h(d) = 2, but the 
non-equivalent forms simply differ by a sign. 

5. The matrices 


A, = 


1 

6 > 

I and A 2 = ( 

' 7 


\ 6 

37 j 


v 12 

31 ; 


correspond to d = 40 with D = 6 2 d, l(Af) = 6. Then Qa\ = — 6(x 2 + 6 xy — y 2 ) and 
Qa* = —6(3x 2 +4:xy—2y 2 ) are the two corresponding (non-trivially) non-equivalent 
forms. 
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Remark. In the case when h(d) is larger than 1 in general we do not have a 
simple formula for the multiplicities for a particular SoZ-manifold but only for the 
disjoint union of SoZ-manifolds with non-equivalent forms of given discriminant d. 

The fact that the multiplicities are large and not sensitive to the change of 
the parameters in the metric seems to be remarkable. A possible explanation of 
the rigidity of multiplicities for SoZ-manifolds is in the hyperbolicity hidden in the 
topology of the manifolds. 

Remark. The same numbers Nq^u) appear in the harmonic analysis on Sol- 
manifolds as the multiplicities of the irreducible Sol -representations in C°°(M\) 
(see Chapter 1 in |S]). Although this fact has a similar origin it does not explain 
the degeneracy of the spectrum of A. In fact one can check that that the same 
degeneracy holds for a more general class of the metrics on 

(22) ds 2 = a(z) dx 2 + 2/3(z)dxdy + ~f(z)dy 2 + dz 2 


where 


( a ( z ) /?(-) 
V /?(-) 7 (a) 


^ = exp(— zB) t 


a 0 (z) /3 0 (z) 
Po{z) 7o(*) 


exp(— zB). 


Here exp B = A and ao(z), (3q{z), 70 ( 2 ) are arbitrary real 1-periodic functions with 
the only condition that the form ds 2 = oto(z)dx 2 + 2/3 0 (z)dxdy + 'yo(z)dy 2 is pos¬ 
itive definite for all 2 . The SoZ-invariant metrics © correspond to the case when 
these functions are constant. The degeneracy in this case follows again from the 
separation of variables and thus is not directly related to the SoZ-invariance. 

It is interesting to compare the Sol -case with the spectra of flat tori T 2 . It is 
easy to see that in the last case the answer will depend drastically on the metric 
parameters (or equivalently, on the geometry of the basic parallelogram). For ex¬ 
ample, if it is a square then the spectrum up to a multiple is given by the values 
of the standard quadratic form n = x 2 + y 2 , and the multiplicity are given by the 
Gauss’ famous formula 


m = 4(IVi(n) - N 3 (n)), 


where Ni(n) and N 3 (n) are the numbers of the divisors of n with the residues 1 
and 3 modulo 4 respectively. If however the basic parallelogram is generic then all 
multiplicities are 2 (which is due to the central symmetry of the problem). 


7. SEMICLASSICAL ANALYSIS AND WEYL’S LAW 


It is instructive to see how our exact calculation of the spectrum agrees with the 
famous Weyl’s law m, which says that for a quantum system the number N(A) of 
the eigenvalues £ < A for large A asymptotically is equal (up to a factor ( 2n )~ n ) 
to the volume of the domain in the classical phase space with the energy less than 
A. For our Laplace-Bcltrami operator m this means that 


(23) 


N( A) 


-ttA W Vd ( M p 
3 (2t r ) 3 


3 (27t) 3 


where Vol(M^) is the volume of our SoZ-manifold (which equals the area of the 
fibre A(T 2 ) since the length in 2 -direction was assumed to be 1). 

Let us count now the eigenvalues £ using the results of section 5. Let us assume 
for simplicity that cos 9 — 0, so besides the trivial part they coincide with the 
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eigenvalues of the Mathieu operator 


(24) 

where as before 


M = - 


d? 2 


\v\ cosh 2 pz 


(25) 


87r 2 Q 

VdA(T 2 ) sing 


and Q is the corresponding binary quadratic form. 

First of all let us use the well-known fact from number theory (see e.g. El) that 
for large Q o the number of lattice points (modulo A) with values of \Q\ less than 
Q o is proportional to the area of the fundamental domain up to Q o : 


M(Q 0 ) ~ 4-^=Q 0 • 

The factor of four counts lattice points related by the symmetry given by chang¬ 
ing the sign of both p u and p v and also accounts for the states in the quadrants 
where Q is of opposite sign. 

For fixed value of Q (hence v) there is a whole line of eigenvalues of the Mathieu 
operator m- The number of these eigenvalues up to energy A for large A is given 
asymptotically by the action integral 

/(A, Q) = — a a/A— \v\ cosh 2gz)dz, 


which is of course the area of the domain in the phase plane with energy less than 
A divided by 2n. This can be simplified to 


27r/r/= \/A^ y/T^gcosh 2 ^d(), g= M/A. 

With £ = cosh(2£) this becomes a standard elliptic integral (see e.g. i501 1 


(26) 


27 ip 

vz 


1 = 2 


rl/9 VT^ 


d£ = 4 y/l + g(K(k) - E(k)), k 2 = 


1 -g 
1 +9 


Let us denote this expression f(g). 

Thus we see that the total number of states up to energy A is 

N{ A)~ J M\Q)I(A,Q)dQ = 2^ -^=f(g)dg, 

o 2 

where C = • s - integral over g is best performed by treating it as 

a double integral over g and £. Introducing rj = gt; and performing the 77 -integral 
first gives 

iW if-s—u. 


3 


Thus we have 


(27) N(A)~A 

which agrees with Weyl’s formula (El when 0 = tt/2. 

For general 6 we have £ = A — vcosd. In that case we need to compute 


X± 


y/1 — (/(cosh 2z ± cos 9)dzdg 
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E 



Figure 3. First 15 states of the cosh-Mathieu equation in depen¬ 
dence of the parameter. 


over the domain 0 < g < 1/(1 ±cos0) and \z\ < z$ where zq is the smallest positive 
root of the integrant. The transformation 

r) = g(cosh2,z + cos 6), £ = cosh 2z 


folds the integration region to the rectangle £ > 1 and 0 < rj < 1 and the integral 
becomes _ 

X±=2 f f - x/1 ~ V - dr/dg . 

J J (£ ± cos 9)\/t; 2 — 1 

The integral over r] is easily done as before while the integral over £ evaluates to 


Y 4f±(0-§) 

A ± — o--' 

3 sm 9 

Hence the sum of the contributions from the two cases of positive and negative Q 
is 


sin0(X + + X_) 


4 

— 7 T 


3 


as before. This computation shows that 9 determines the relative number of states 
between the regions with positive and negative Q , namely X + /X_ = 9/(tt — 8). 

Let us look what this calculation gives for the first eigenvalues of A. There are 
two opposite cases depending on whether geometric parameter A = A(T 2 ) sin 9 is 
small or large. Let us assume again for simplicity that 9 = n/2, then A is simply 
the area of the fibre. 

The small A corresponds to the ”rope-like” ^oZ-manifolds. In this case the first 
eigenvalues are ’’trivial”: Ek = Ak 2 n 2 , k = 0,1, 2... which correspond to Q = 0. 

The second case when A is large is more interesting. In that case the parameter 
v in the Mathieu operator is small. The action integral for small v has the 
asymptotics 1 


I 


VX A_ 

HTT 11 H' 


This suggests the following asymptotics for the eigenvalues for small v: 

(fink ) 2 


A k = 


(In M) : 


r, * = 0,1,2.... 


^The fact that the action diverges logarithmically for Q —> 0 seems to be surprising. To our 
knowledge this is the first example of Liouville integrable system on a compact manifold for which 
the action diverges on approach of a singular level (with energy fixed). 
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Note that although Afc —> 0 as v —> 0 the decay is slower than any power of v. The 
first 15 states are shown in fig. 0 

For the corresponding Sol-manifolds this gives the following behaviour of the 
first eigenvalues 

(28) 4 ,[ 7 ] ~ ( ln |(7Q A »([ 7 ])|)2’ * = 0,1,2,..., b]er*\{0}/T 

O 2 

for small C = — ,Th --. If we order these eigenvalues £q = 0 < £\ < £ 2 , ■■ ■ 

\/DA{T Z ) sin 6 

then we have 


(29) 


£, 


(f™ ) 2 ( , O ln Q) S 
(InC ) 21 In C J ’ 


j = 1,2 ,..., 


where Qo = 0 < Qi < Q 2 , ■ • • are positive values of the form Qa * listed in increasing 
order and we have assumed that Qj <C In particular, 


(/it) 2 £j - £ x ln(Qj/Qi) 

(ln C) 2 ’ £< - £i ~ MQiM) ’ 


so when C is small we ’’see” the values of the quadratic form Qa * straight from 
the spectrum. 

Note that the question about the next order term in Weyl’s law is non-trivial. 
For the simpler case of iViTmanifolds some results in this direction can be found 

in 20]. 

One can also look at the corresponding Minakshisundaram-Plejel asymptotic ex¬ 
pansion, which is an important characteristics of the spectra (see e.g. D3H)- In par¬ 
ticular the second coefficient in this expansion is proportional to the integral of the 
scalar curvature K. A straightforward calculation shows that for the SoZ-manifold 
M\ the principal sectional curvatures are ± sin 26 log 2 A and —2 (sinolog A) 2 , so 


K = —2(sin 6 log A) 2 


and thus is always negative. 


8. Spectral Statistics 

The spectral statistics of integrable and chaotic systems is quite different, see, 
eg., mmm- As we have seen the geodesic flow on Sol-manifolds has properties 
of both, integrable and chaotic systems. Therefore it is a natural question what 
the spectral statistics of the Sol-manifolds is like. Note that according to the 
Berry-Tabor conjecture [2J integrable systems should have Poisson distributed level 
spacing. We are going to show that this is not the case for Sol-manifolds. 

The reason is the high multiplicities of the eigenvalues. Indeed, for the simplest 
positive quadratic form Q 0 = x 2 + y 2 the classical result due to E. Landau says 
that the number of integers up to a number K represented by this form grows as 
K /Vlog K. If there would be no degeneracies then this number would grow like 
the area of the fundamental region, which is proportional to K. This means that 
most of the level spacings of the values of Q 0 are zero. 

According to P. Sarnak EH El. a similar fact is true for indefinite forms as well, 
namely 

The number of positive integers up to K that can be represented by a given indefinite 
quadratic form Q grows not faster than 0{K/y/\og K). 
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Qmax = 10 A 5., # = 21749 % Qmax = 10 A 5., # = 21749 




Qmax = 10*7., # = 2154323 


Qmax = 10*7., # = 2154323 



0.1 



Figure 4. Level spacing statistics of the indefinite binary qua¬ 
dratic form Q(x,y) = —x 2 + xy + y 2 (left) with degeneracies re¬ 
moved (right) 


Combining this with the results of section 5 we have the following 

Theorem 4. The level spacing distribution for the spectrum of Sol-manifolds M\ 
is not Poisson and hence the Berry-Tabor conjecture does not hold in this case. 

This is particularly interesting because this statement is not sensitive to change 
of the metric in the Sol- class © ( or even more general class (J22J0- 

Let us illustrate this in the example of the cat-map A. In fig. 2] (left) the level 
spacing statistics for the indefinite binary quadratic form Qa = ~x 2 + xy + y 2 is 
shown for three different values of Q max . 

Since the cat-map is the product of two involutions, there is a simple reflection 
symmetry in the lattice, which causes almost all states to be at least twofold de¬ 
generate. This discrete symmetry needs to be factored out before the level spacing 
statistics can be studied. The involutions Ri with R 2 = Id are 

a=r 2 r 1 , ri=(\ _ y , :;). 

The fixed line of R± is the line y = —a;/2, and factoring the fundamental region 
in fig. m by R\ simply cuts the fundamental region in half along this line. Since 
the values of Q are integers we chose to present the raw level spacing statistics, 
i.e. without unfolding the spectrum first. The number of lattice points found up to 
the corresponding Qmax in the reduced fundamental region is given in the heading 
of each figure, and the ratio approaches lnA/(2 y/d). The figures clearly show that 
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Figure 5. Image of the lattice 1? in ( p x ,p y ) under the momentum 
map (i*i, F 2 ) for fixed energy where A is the cat-map. Left: origin 
at the centre. Right: distorted standard lattice away from the 
origin. 


the proportion of degenerate levels grows in agreement with what we said above. 
When the degenerate levels are discarded in the statistics fig. 01 (right) shows that 
the distribution appears to converge to some non-universal shape. We would like to 
mention here the paper m where the moments of the intervals between the sums 
of two squares were studied. 

9. Quantum monodromy 

In view of the previous results the appearance of quantum monodromy in our 
problem is quite natural. However there is a problem with this notion in our case 
which we want to discuss first. 

As it was shown in J2J the geodesic flow on Sol-manifolds can not have three 
analytic integrals (see Section 3 above). A similar fact holds in the quantum case. 
Namely, one can show that the algebra of the differential operators on the Sol- 
manifold M\ commuting with the Laplace-Beltrami operator A 111 511 is generated 
by A and This means that our quantum problem does not have enough 

quantum integrals, at least in the class of the differential operators and therefore 
it is not clear if we can apply the rigorous treatment of quantum monodromy from 

HU¬ 
SO in this section we will treat the quantum monodromy on Sol-manifolds on 
the intuitive level paying more attention to geometry rather than to analysis. As 
we have already shown, the set of eigenfunctions is in a natural one-to-one corre¬ 
spondence with T*/A* x N, [ 7 ] £ T*/A*, k £ N. The fundamental domain of A* 
is shown in fig. [21 It is natural to represent the orbit space as a lattice on the 
cone obtained by gluing the edges of the fundamental domain of A* on the plane 
(more precisely we should consider four different cones corresponding exactly to 
four families of Liouville tori). 

Quantum monodromy arises when we pass around the vertex of the cone. It 
is clear that the basis of the lattice will undergo the transformation A. On the 
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Figure 6. Image of the lattice 1? in ( p x ,p y ) under the momentum 
map (Fi, F 2 ) for fixed energy and A the cat-map. Left: Sol-flower 
with \Q\ < 60 2 . Right: Sol-flower with \Q\ < 90 2 


other hand, nothing happens to the third direction corresponding to the parameter 
k. Therefore the quantum monodronry for the Sol-manifold M\ is given by the 
matrix 

A* 0 
0 1 

We want to emphasize that in this case quantum monodromy has a purely topo¬ 
logical nature. It is determined by the topology of the underlying manifold, and 
not by properties of the metric. It does not depend on the parameters E , G, F, 
moreover the monodromy remains the same for all metrics of the form 

ds 2 = ds 2 + dz 2 

where ds 2 is a flat metric on fibres T 2 with coefficients depending on z. In the 
previously known examples (like the geodesic flow on the 3-dimensional ellipsoid of 
revolution, |28l l the metric g played the principal role. 

In figures 13 El and Q we demonstrate the quantum monodromy of the Sol- 
manifold M\ related to the cat-map 


A = 




To make the image of the lattice uniform we have slightly modified the classical 
integrals f \, f 2 from section 3 as follows: 

F\ = \/\Q\ cos 27r A 

F 2 = VW\ sin 

where Q = p u p v and f3 = ' . The image of the lattice under the map F we 

call Sol-flower (see fig. El Note that an alternative choice (3 = would give 
a similar picture and the freedom of the rescaling of the eigenvectors e u ,e v leads 
simply to a rotation of the plane (Fi,F 2 ). Fig. El illustrates that away from the 
origin the lattice is simply a deformed standard lattice. 
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Figure 7. Parallel transport of basic directions in the image of 
the momentum map. 


A nice property of the map (p u ,Pv) —> (Fi, Fz) is that it changes the area simply 
by a constant multiple: it is easy to check that 

dFi A d F 2 = 7^--dpu A d p v . 

In A 

It is interesting to mention that the multiplicity problem becomes the standard 
’’circle problem” if one replaces the square lattice by the 5oZ-flower (but of course 
it does not help to compute them). 

When a fundamental cell is chosen in the FoZ-flower as indicated in grey in 
fig-0 the monodromy can be observed as follows: A line extending a basis vector 
is parallel transported in the lattice. After completing a cycle about the origin 
this direction is changed. The left picture shows images of the lines (p x ,p y ) = 
(30 — 21, —j — /), j = 0..27, l = 0..5. The right picture shows images of the lines 
( Px,P y ) = (30 — l,—j — l), j = 0..30, l = 0..5. Denote the direction of the line 
shown in the left part of fig. 0 by ei, and the one on the right part by e 2 - The 
preimages of these basis vectors in fig. dare — {2e x + e y ) and — (e x + e y ). Parallel 
transporting e\ clockwise by increasing j gives e\ + e 2 (determining the second 
row of A), while parallel transporting e 2 counterclockwise by decreasing j gives 
—ei + e 2 (determining the first row of A -1 ). Since A £ SL( 2,Z) this determines 
the cat-map. 


10. Concluding remarks 

The SoZ-geometry from a dynamical point of view has the special property of 
being on the border between integrability and chaos. Integrability is reflected in 
the solvability of the corresponding group while the chaos is related to a hidden 
(partial) hyperbolicity. This makes the Sol -case of particular interest and explains 
why the geodesic problem on the SoZ-manifolds has both integrable and chaotic 
features. As we have seen the quantum case gives a new interesting twist to the 
story by bringing arithmetic into play. 

Atiyah, Donelly and Singer [Q considered a more general case of SoZ-manifolds 
which are T n+l torus fibres over T n . Much of our analysis can be generalised to 





26 


A.V. BOLSINOV, H.R. DULLIN, AND A.P. VESELOV 



Figure 8. Slice of 3?($[ 7 ] ) *) at x = 0 for 7 = (1,0), k = 15 (left) 
and 7 = (12, —5), k = 2 (right) for the cat-map. The appearance 
of the eigenfunction whose Q( 7 ) is small and thus close to the 
classical chaos is more irregular. 


this case as well. The quantum Toda lattice Hamiltonian will appear then as a 
generalisation of the modified Mathieu operator. Some very interesting results in 
the corresponding classical problem were found recently by Leo Butler in [E| . 

It would be also interesting to study in more detail how the chaos (at the degener¬ 
ate level Q = 0) of the classical system manifests itself in the quantum version. We 
showed that the spectral statistics provides a counterexample to the Berry-Tabor 
conjecture, but it cannot be taken as an indicator of chaos. One simple observation 
is that the trivial eigenfunctions <f>o lS are asymptotically ‘uniformly distributed’ 
on the manifold. Hence the subset of eigenfunctions that are associated with the 
classical chaos are quantum unique ergodic, cf. EH Already at relatively small 
quantum numbers it can be seen that the nodal lines are more complicated when 
Q is small, see fig. |H1 
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